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The normal state in the hole underdoped copper oxide superconductors has proven 
to be a source of mystery for decades. The measurement of a small Fermi surface 
by quantum oscillations on suppression of superconductivity by high applied magnetic 
fields, together with complementary spectroscopic measurements in the hole under¬ 
doped copper oxide superconductors, point to a nodal electron pocket from charge 
order in YBa2Cu3 06 +a. Here we report quantum oscillation measurements in the 
closely related stoichiometric material YBa 2 Cu 408 , which reveal similar Fermi surface 
properties to YBa2Cu3 06+5, despite an absence of charge order signatures in the same 
spectroscopic techniques such as x-ray diffraction that revealed signatures of charge or¬ 
der in YBa2Cu3 06 + 4 . Fermi surface reconstruction in YBa2Cu4 08 is suggested to occur 
from magnetic field enhancement of charge order that is rendered fragile in zero mag¬ 
netic fields because of its potential unconventional symmetry, and/or its occurrence as 
a subsidiary to more robust underlying electronic correlations. 

PACS numbers: 71.45.Lr, 71.20.Ps, 71.18.+y 


Significance 

Results from quantum oscillation and spectro¬ 
scopic measurements have suggested charge order 
as responsible for the creation of a nodal electron 
pocket in the YBa2Cu306+A and HgBa2Cu04+a 
families of underdoped cuprates. However, the 
situation in the pristine YBa2Cu408 family remains 
ambiguous. Our high-precision quantum oscillation 
measurements point to a very similar nodal elec¬ 
tron pocket in the stoichiometric cuprate family 
YBa2Cu408, despite the non-observation of charge 
order in diffraction experiments. Our findings in 
YBa2Cu408 indicate Fermi surface reconstruction 
in the hole underdoped cuprates associated with 
the magnetic field enhancement of charge order, 
the fragile low-magnetic field character of which is 
reflected in the short correlation length reported in 
various materials families. 


Introduction 

The normal state of the underdoped copper oxide su¬ 
perconductors has proven to be even more perplexing 
than the d-wave superconducting state in these materials. 
At high temperatures in zero magnetic fields, the normal 
state of the underdoped cuprates comprises an unconven¬ 
tional Fermi surface of truncated ‘Fermi arcs’ in momen¬ 


tum space, which is referred to as the pseudogap state [1]. 
At low temperatures in high magnetic fields, quantum os¬ 
cillations reveal the non-superconducting ground state in 
various families of underdoped hole-doped copper-oxide 
superconductors to comprise small Fermi surface pock¬ 
ets [2-15]. These small Fermi pockets in YBa2Cu306+6 
have been identified as nodal electron pockets [2, 3, 
11, 16, 17] originating from Fermi surface reconstruc¬ 
tion associated with charge order measured by x-ray 
diffraction [18-20], ultrasound [21], nuclear magnetic res¬ 
onance [22], and optical reflectometry [23]. However, var¬ 
ious aspects of the underlying charge order and the as¬ 
sociated Fermi surface reconstruction remain obscure. A 
central question pertains to the origin of this charge or¬ 
der, curious features of which include a short correlation 
length in zero magnetic held that grows with increas¬ 
ing magnetic field and decreasing temperature [20]. It 
is crucial to understand the nature of this ground state 
order that is related to the high temperature pseudogap 
state and delicately balanced with the superconducting 
ground state. Here, we shed light on the nature of this 
state by performing extended magnetic held, tempera¬ 
ture, and tilt angle-resolved quantum oscillation experi¬ 
ments in the stoichiometric copper oxide superconductor 
YBa2Cu408 [24]. This material with double CuO chains 
has fixed oxygen stoichiometry, making it a model system 
to study. YBa2Cu408 avoids disorder associated with 
the fractional oxygen stoichiometry in the YBa2Cu306+5 
chains, which has been shown by microwave conductivity 
to be the dominant source of weak-limit (Born) scatter¬ 
ing [25]. 

Intriguingly, we find magnetic held- and angle- 
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dependent signatures of quantum oscillations in 
YBa2Cu408 [13, 14] that are very similar to those in 
YBa2Cu306+i, indicating a similar nodal Fermi surface 
that arises from Fermi surface reconstruction by charge 
order with orthogonal wave vectors [16]. However, the 
same x-ray diffraction measurements that show a Bragg 
peak characteristic of charge order in YBa 2 Cu 306 +i for a 
range of hole dopings from 0.084 < p < 0.164 [19, 20, 26[ 
hoave, thus far, not revealed a Bragg peak in the case of 
YBa2Cu408 (Fig. 3 in ref. [19]). We suggest that charge 
order enhanced by applied magnetic fields reconstructs 
the Fermi surface in YBa2Cu408, whereas charge order 
is revealed even in zero magnetic fields in YBa 2 Cu 306 +i 
because of pinning by increased disorder from oxygen 
vacancies. 


Results 


a c 



Fig. 1 shows quantum oscillations in contactless con¬ 
ductivity [27] measured up to 90 T in YBa2Cu408 and at 
different temperatures from 1.3 K to 8.0 K. The extended 
magnetic field range and increased sensitivity compared 
with previous quantum oscillation experiments [13, 14[ 
enable precision measurements of the quantum oscilla¬ 
tion frequency spectrum and effective quasiparticle mass 
of YBa2Cu408. 

Preliminary quantum oscillation measurements on 
YBa2Cu408 accessed two [13] to four [14] oscillation pe¬ 
riods over a restricted magnetic field range, for a mag¬ 
netic field angle parallel to the crystalline c-axis. A sin¬ 
gle quantum oscillation frequency of 660 ± 30 T was 
reported, whereas the scatter of the quantum oscilla¬ 
tion amplitude as a function of temperature precluded 
a determination as to whether a Lifshitz-Kosevich form 
is obeyed, or an accurate extraction of a quasiparticle 
effective mass. Our present quantum oscillation mea¬ 
surements over an extended magnetic field range access 
more than seven oscillation periods, revealing for the 
first time a pronounced quantum oscillation beat struc¬ 
ture characteristic of multiple frequencies, very similar 
to YBa 2 Cu 3 06.56 [6, 9, 28], with a dominant frequency 
of 640 T. Our precision measurements of quantum os¬ 
cillation amplitude as a function of temperature shown 
in Fig. lb and c further reveal a distinctive Lifshitz- 
Kosevich form, characteristic of Fermi Dirac statistics. 
A fit to the Lifshitz-Kosevich form yields a quasiparticle 
effective mass of 1.8 ± IrUe in YBa 2 Cu 408 , which is, in 
fact, very similar to that measured for YBa 2 Cu 3 06 , 56 . 

Angular measurements as a function of tilt angle to 
the applied magnetic field are required in order to iden¬ 
tify the origin of the multiple frequency spectrum that 
we observe in YBa2Cu408. Fig. 2 shows the quantum 
oscillations in YBa 2 Cu 408 measured up to a maximum 
tilt angle of 0 « 56°. A few key features are notable. 
First, the beat pattern at the zero degree tilt angle per¬ 


FIG. 1: (a) Quantum oscillations measured in the contact¬ 
less resistivity of YBa 2 Cu 408 for the magnetic field parallel 
to c-axis. (b) Fourier transform of the measured quantum 
oscillations at different temperatures, (c) Plot of quantum 
oscillation amplitude as a function of temperature (symbols) 
accompanied by a Lifshitz-Kosevich fit (red line). 


sists up to high tilt angles. Second, the prominent Ya- 
maji amplitude resonance [29] expected at 6 « 52° for 
a neck-belly warped Fermi surface in YBa2Cu408 is ab¬ 
sent up to the measured high tilt angles. The absence of a 
prominent Yamaji angle in YBa2Cu408 enables us to rule 
out a Fermi surface where the observed frequency spread 
arises from a dominant neck and belly fundamental warp¬ 
ing. Instead, both of the observed features are consistent 
with a quasi-two dimensional Fermi surface in which the 
multiple frequency spectrum originates from a splitting 
of frequencies as opposed to a fundamental neck and belly 
warping. Such a quasi-two dimensional Fermi surface is 
similar to that identified in YBa 2 Cu 306 +^ [16, 28[. 

Another clue as to the origin of the observed quan¬ 
tum oscillations in YBa 2 Cu 408 is obtained by inspecting 
the evolution of the amplitude and phase of the quan¬ 
tum oscillations as a function of tilt angle. Fig. 3b 
shows the cross-correlation function between the mea¬ 
sured quantum oscillations in YBa 2 Cu 408 and a phase- 
matched sinusoidal function with frequency 640 T, which 
is averaged over the indicated magnetic field range, re¬ 
ferred to as the correlator. We find that the corre¬ 
lator is very similar to that previously measured for 
YBa 2 Cu 3 06.56 [11, 16] and shown in Fig. 3a. The corre¬ 
lator for YBa 2 Cu 408 reveals a zero crossing and phase 
inversion of the quantum oscillation amplitude, signalling 
a spin zero at a tilt angle of 9 Ri48°. A spin zero feature 
arises on destructive interference between quantum oscil¬ 
lations from two spin channels at certain special angles. 
At these spin zero angles, the ratio of the Zeeman energy 
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FIG. 2: Quantum oscillations (in colour) measured in the 
contactless resistivity of YBa 2 Cu 40 g for different angles of 
inclination ( 6 ) of the magnetic field to the c-axis plotted as a 
function of Simulated quantum oscillations (in gray) 

of a bilayer split quasi-two dimensional Fermi surface model 
shown in Fig. 4 [16] are shown for the parameters listed in 
Table I. The inset shows a schematic of the crystal tilt angle 
to the magnetic field. 


to the cyclotron energy — captured by the spin splitting 
factor = cos „ 1 — crosses zero, and inverts 

sign. Spin zero angles are located at 


dsz = cos 


(2n -|- l)m,e 


( 1 ) 


where n is an integer, g* is the effective p-factor, assumed 
here to be isotropic for simplicity (see Methods and Ta¬ 
ble I for the more general case), rrie is the electron mass, 
and mj| is the effective quasiparticle mass for B parallel 
to the crystalline c-axis [30]. 

A value of rn*^g* can be extracted by the measurement 
of at least two spin zeros, given the two unknown quan¬ 
tities n and g*, as was done for YBa 2 Cu 3 06.56 [16]. To 
interpret the single spin zero observed at « 48° 

in YBa 2 Cu 40 g, a comparison with YBa 2 Cu 3 06.56 is in¬ 
structive, in which case two spin zeros were observed at 
^first ~ 55° [10, 16[ and 9^1^ « 66° [16]. Good agree¬ 
ment with experimental data is yielded by associating 


an index value of n = 2 with the first spin zero, for pa¬ 
rameters including anisotropic ^-factors (see Methods) 
given in Table I. Parameters obtained for YBa2Cu408 
are remarkably similar to those obtained in the case of 
YBa 2 Cu 3 06.56 [16], and are in good agreement with ex¬ 
perimental data (solid lines in Fig. 3). 


a 
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FIG. 3: a Symbols represent the cross-correlation function be¬ 
tween the measured quantum oscillations in YBa 2 Cu 306.56 
and a phase-matched sinusoidal function with frequency 
534 T, which is averaged over the indicated magnetic field 
range (data from ref. [16]). b Symbols represent the cross¬ 
correlation function between the measured quantum oscilla¬ 
tions in YBa 2 Cu 408 and a phase-matched sinusoidal function 
with frequency 640 T, which is averaged over the indicated 
experimental magnetic field range. Solid lines show the fits 
using the equation and parameters in Table I. 

Given the similarities we find between YBa2Cu408 
and YBa 2 Gu 306+6 in terms of the measured multi¬ 
ple quantum oscillation frequency spectrum, split quasi- 
two dimensional Fermi surface, and spin zero angles 
revealed by the correlator, we compare the measured 
quantum oscillations in YBa2Gu408 with the Fermi sur¬ 
face model fit to the measured quantum oscillations in 
YBa2Cu306-i-i [11, 16, 17]. We consider a model in 
which cyclotron orbits are associated with a nodal bi¬ 
layer split Fermi surface from charge ordering shown in 
Fig. 4. Here the quantum oscillation frequency spread 
is associated with a splitting of the Fermi surface aris¬ 
ing from tunnelling between bilayers. Magnetic break¬ 
down tunnelling in the nodal region where the splitting 
is smallest can then give rise to a series of combination 
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Parameter 

Description 

YB32CU408 

YBa2CU306.56 

Fo 

Quantum oscillation frequency 

639 T 

534 T 

^ ^twofold 

Staggered twofold warping frequency 

— 

15 T 

split 

Bilayer splitting frequency 

91 T 

90 T 

m|, 

Quasiparticle effective mass 

1.8 rPe (fixed) 

1.6 rrig (fixed) 

So 

Magnetic breakdown field 

4.2 T 

2.7 T 

9|d 

g-Factor 1 

2.0 

2.1 

9110 

g-Factor 2 

0.1 

0.4 


g-Factor anisotropy 1 

1.6 

1.4 

fo 

g-Factor anisotropy 2 

0.8 

0.2 

4> 

Phase 

-1.6" 

0 (fixed) 


TABLE I: Parameters used to simulate the oscillatory waveform for a quasi-two dimensional split Fermi surface model shown in 

Fig. 4, represented by the Eqn. 4 / twofold ~ NjIRmbRsRdRt]] cos ^ Jcosg + • Here, iiMB is the magnetic breakdown 

amplitude reduction factor (defined in Methods), and Nj counts the number of instances that the same orbit is repeated within 
the magnetic breakdown network. Ru is the Dingle damping factor, Rt is the thermal damping factor, and Rg is the spin 
damping factor (defined in Methods). The tabulated values used to simulate the quantum oscillation waveform yield good 
agreement with experiment as a function of B and 8 (Figs. 2 and 3). The effective mass is taken to be a fixed quantity, having 
been determined independently from temperature-dependent measurements [9]. The parameters are the same for all of the 
orbits, except for those denoted by subscripts □ and <(>, which each correspond to a subset of orbits as defined in the text. The 
values of and here represent parameters used for the simulation rather than unique identifications. The parameters used 
for YBa 2 Cu 306.56 shown for comparison are taken from ref. [16]. 


frequencies, which is discussed in refs. [16, 28]. Figs. 2 
and 3 shows that this model can simulate the angular and 
magnetic field dependence of the quantum oscillations 
measured in YBa2Cu408 reasonably well. The parame¬ 
ters used for the model simulation are in fact very similar 
to those used to simulate quantum oscillations measured 
in YBa 2 Cu 3 06.56 (Table I) [16]. The small staggered 
twofold warping included in the case of YBa 2 Cu 3 06.56 is 
not included for YBa 2 Cu 408 , given the restricted angu¬ 
lar range over which quantum oscillations can be accessed 
and the likely weaker amplitude warping associated with 
a longer c-axis. We note that there is some deviation 
from the model at a few of the highest angles. This devi¬ 
ation may be caused by additional lifting of degeneracy 
of frequency components from effects such as a distortion 
of the simple tetragonal crystal structure or other details 
of Fermi surface geometry beyond those considered in the 
present model. 


Discussion 

Our findings therefore reveal closely similar quan¬ 
tum oscillation features in YBa 2 Cu 408 compared to 
YBa 2 Cu 3 06 , 56 , showing (i) a similar quasiparticle effec¬ 
tive mass from a Lifshitz-Kosevich fit to the amplitude 
dependence as a function of temperature, (ii) a spread of 
multiple frequencies yielding a prominent beat structure, 
(iii) an angular dependence of quantum oscillation fre¬ 
quencies consistent with Fermi surface splitting accom¬ 
panied by magnetic breakdown rather than fundamen- 





FIG. 4: Brillouin zone cross-section, showing a schematic of a 
nodal electron pocket created by charge ordering wave vectors 
Qi = 27r(±^,0) and Qa = 27r(0,±^) [11, 16, 17]. a shows 
the two Fermi surface cross-sections of frequency Fi = Fq — 
2AFspiit and Fe = Fq-I- 2AFspia. The gap separating bonding 
and antibonding surfaces is expected to be smallest at the 
nodes [31]. A cut through the kz = 0 plane of the Brillouin 
zone shows the possible magnetic breakdown orbits [16, 28]. 
b — e show the range of possible magnetic breakdown orbits 
(Fa = Fo — AFspiit, F 3 = Fq, F 4 = Fq and Fg = Fq -|- AFspia) 
as listed in Table 1. 
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tal neck and belly warping, and (iv) a similar valne of 
m*g* from the angnlar dependence and spin zero ob¬ 
served in the correlator. The common Fermi surface fea- 
tnres point to the same origin of Fermi surface recon¬ 
struction in YBa 2 Cu 408 as in YBa 2 Cn 306 + 5 , which is 
particularly relevant given the non-observation of charge 
order signatures thus far in YBa2Cn408. In the case of 
YBa 2 Cn 306 + 5 , indications are that the electronic struc- 
tnre comprises a nodal electron pocket from Fermi sur¬ 
face reconstruction by a superstructure with orthogonal 
wave vectors Qi = 27r(±^,0) and Q 2 = 27r(0,±^) as¬ 
sociated with charge order observed by x-ray diffraction 
and other experiments (where « J 2 ~ 0.3) [19, 20, 22]. 

An unusual aspect of the charge order measured in 
YBa 2 Cu 306+5 is the short average correlation length, 
which has been measured to be of the order of 100 A 
in the basal plane and the order of 10 A along the c-axis 
from x-ray diffraction measurements in a magnetic field 
of 17 T [20]. The measured average correlation length in 
the basal plane is comparable with the cyclotron radius 
of 100 — 200 A in the magnetic field range of 30 — 100 T 
where quantum oscillations are measured [2-12]. More¬ 
over, given that the anisotropic Fermi velocity associated 
with the observed pocket has a magnitude along the c- 
axis that is approximately two orders of magnitude lower 
than in the basal plane [28, 31, 32], the component of the 
cyclotron trajectors along the c-axis is also well within 
the c-axis correlation length. The measured charge or¬ 
dering wave vectors with finite correlation length inferred 
from x-ray diffraction experiments represents an average 
over all parts of the sample, thus representing a lower 
bound on the correlation length associate with the most 
highly ordered regions of the sample, which are accessed, 
for instance, by quantum oscillation measurements. It 
has remained unclear thus far whether the reported short 
correlation length is due to the nucleation and pinning of 
charge order at impurity/defect sites, or whether it is be¬ 
cause of the disruptive effects of impurities/defects in an 
intrinsically long range charge ordered state ]33, 34]. 

Our findings in YBa2Cu408, where impurity pinning 
centres are reduced compared to YBa 2 Cu 306 +^ because 
of the fixed oxygen stoichiometry, suggest that charge 
order of inherently short correlation length arises in zero 
magnetic fields in YBa 2 Cu 306+5 on account of impurity 
sites that act as pinning potential centres. Charge order 
in the well-known charge density wave material NbSe 2 
has, for instance, been shown to initially develop around 
impurities ]35]. In the case of YBa2Cu408, however, 
where such impurity pinning centres are reduced, applied 
magnetic fields are required to tilt the balance of energy 
scales such that charge susceptibility is further enhanced, 
and charge order is revealed to reconstruct the Fermi sur¬ 
face. 

Our findings point to a fragile form of charge order 
in the non-superconducting ground state of the under¬ 
doped cuprates which is enhanced by applied magnetic 


fields to yield Fermi surface reconstruction. An inter¬ 
esting question pertains to the origin of such intrinsi¬ 
cally fragile charge order. A contributing factor may be 
the potential unconventional character of the observed 
charge correlations ]36]. Furthermore, charge order may 
arise as a corollary to more robust correlations in the un¬ 
derdoped cuprates, such as underlying strong spin cor¬ 
relations [37]. The interplay between these correlations 
may be manifested as other order parameters such as 
Amperean order ]38], pair density wave order ]39], d- 
wave checkerboard order ]40, 41], and quadrupolar or¬ 
der [42, 43], which have also been proposed to appear in 
fluctuating form above the superconducting temperature. 

Methods 

Quantum oscillation simulations include conventional 
thermal, Dingle, and spin damping factors of the same 
form used for previous comparisons with quantum oscil¬ 
lations measured in the underdoped cuprates and other 
layered families of materials [30]. 

The thermal damping factor is given by 


sinh Xj 

(where Xj = 2'K‘^kBmg^T/heB), is the Boltzmann fac¬ 
tor, and T is the temperature; mgj = ^|j/ cos 6 (taken to 
be the same for all orbits and the subscript ‘j’ dropped) 
is determined by the projection of B perpendicular to 
the planes (i.e. the projection parallel to the c-axis in 
YBa 2 Cu 306 +a:), and refers to the value of rrigj when 
B is parallel to the crystalline c-axis [taken to be a fixed 
quantity, determined independently from temperature- 
dependence measurements [9]]. 

The Dingle damping factor is given by 

Rd = exp( - Y 

\ B cos 0 J 

where Aj is a damping factor, taken to be the same for 
all orbits, enabling us to drop the subscript j [30]. 

The spin damping factor is given by 



where the anisotropic effective ^-factor has the form ggj = 
g^ y^cos^ 9 + j- sin^ 6. Here, g*^^ refers to the value of gg. 
when B is parallel to the crystalline c-axis, whereas = 
is the anisotropy in the spin susceptibility. Two 
sets of anistropic ^-factors are considered for subsets of 
orbits defined below. Because of the multiple frequencies 
in the model and the restricted angular range measured, 
it is not possible to uniquely identify the g-factors. The 



values of and here represent parameters used for 
the simulation. 

A splitting of the Fermi surface arising from tunnelling 
between bilayers leads to two starting frequencies that 
are denoted as Fi = Fq — 2 i<’spiit and Fg = -Fb + 2 Fspiit. 
Magnetic breakdown tunnelling (in the nodal region 
where the splitting is smallest) gives rise to a series of 
combination frequencies F 2 , F^, F 4 , F^, as discussed in 
refs. [16, 28]. Only two sets of anisotropic p-factors are 
considered; orbits Fi, F 2 , F 4 , F 5 and Fq, which un¬ 
dergo both magnetic breakdown tunnelling and reflec¬ 
tion, are approximated to have the same g-factor g|^ with 
anisotropy ^o, whereas orbit F 3 , which show only mag¬ 
netic breakdown tunnelling without reflection is approx¬ 
imated to have a common ^-factor gy^ with anisotropy 
[16]. The magnetic breakdown amplitude reduction 
factor is given by 

7?mb = 

in which F and F count the number of magnetic break¬ 
down tunnelling and reflection events en route around the 
orbit, having transmitted amplitudes i\/P and — P 
respectively. The magnetic breakdown probability is 
given by P = exp(—Po/F3cos 9), where Bq is the charac¬ 
teristic magnetic breakdown held [30]. 
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